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I. INTRODUCTION 

In the last decade, after the work of Randall and 
Sundrum[l|, |2j, the study of brane- worlds have increased 
considerably, due to its application in explaining funda- 
mentals problems in both particle physics and gravity. 

Besides the early works [H-Q , the importance of brane- 
worlds became apparent almost 20-years later, when the 
hierarchy problem was considered in this context [l| , and 
the same authors realized that a new kind of compact- 
ification was possible Bl , different to that proposed by 
Kaluza and Klein 0, and its generalizations Q . 

On the other hand, branes are the fundamental 
brick for the AdS/CFT correspondence, proposed by 
MaldacenaQ (see also [T(J), which allows to relate the 
spectrum of a string theory in a AdS background with 
the one of a conformal field theory at the boundary of the 
AdS spacetime. Throughout this article only the branc- 
world aspect will be considered. 

Several generalizations of the Randall-Sundrum brane 
worlds have been considered, such as multi-brane 
worlds thick branes or domain walls fl^. more com- 
plicated domain walls solutions with or without a well de- 
fined brane- limit fl3T - [l5| (for an illustration on the brane 
limit see fl6j). localization of gravity in different brane- 



19] 
21] 



world backgrounds 17|, localization of fermions[T 
and their spectrum (20j|. localization of gauge fields 
and many others. 

The aim of this paper is to restrict the moduli space to 
specific variety of domain walls, known in the literature 
as static double walls (33|, found in 13 1 and generalized 
in [3, [l5|. Despite the existence of these generaliza- 
tions, we concentrate ourselves on the original solutions. 
The motivation for studying these static domain walls is 
that they are the minimal parametric generalization of 
the thick version of Randall-Sundrum brane- world [l2j], 
which are suitable phenomenological candidates for al- 
lowing fermion localization. 



In the following section a detailed analysis of the lo- 
calization of fermions is done, for static solutions. Al- 
though the procedure is known and reviewed in several 
articles, new considerations are made and special cases 
considered. In section IIIII the analysis of the spectrum 
and moduli space is done for a family of static double 
domain walls. At the end an appendix is included to 
summarize notation and conventions. 



II. FERMION LOCALIZATION ON STATIC 
D-BRANES 

In the following, thick-branes (or gravitational domain 
walls) are considered. The formal relation between the 
two might be established through a limit [llij . 

This Domain Walls are generated by a scalar field, 
</>(£), which depends on the extra dimension, £ (see j2r|). 
The spacetime has a metric whose line element is given 

by m, 

ds 2 {g) = e 2A ^ (-dt 2 + e 2B ^dx 2 ) + e 2C («>d£ 2 . (1) 



Several Domain Wall solutions are listed on [13Ml5l | . 

The 5-dimensional action for a massive fermion on a 
curved spacetime is given by, 



S 



D 



d 4 xd^y/\g\^ (- X 7 - M) f , 



(2) 



where is the 5-dimensional fermion, the bar indicates 
Dirac conjugation, V is the covariant derivative on the 
curved spacetime, M is the 5-dimensional mass of ^ and 
\g\ is the absolute value of the metric determinant. The 
use of tetrad formalism is implicit, although the notation 
does not suggest it. 

However, it is not possible to localize fermions on this 
set up 22]. A bypass to the problem consist in adding 
an interaction. A Yukawa-like interaction between the 
fermion an the scalar field which generates the domain 
wall does the work. Henceforth, the action to consider 
is, 
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with V is, in general, a polynomial function on cj>. The 
simplest model is V(4>) = (j), and it was considered in 
[13] with a different kind of Domain Wall than those pre- 
sented here. 

In the literature, it has become customary to take into 
account just the case when the 5-dimensional mass of 
the fermion vanishes, i.e. M = 0. However, as will be 
discussed, M could be a rather important parameter for 
selecting plausible phenomenological models. Therefore, 
the distinction between massive and massless fermions is 
made. 



A. Massless 5D Fermions 

In the following section a presentation of the differ- 
ent possible conditions which could be imposed on the 
fermions are give n, as well as the equations of motion of 
their profiles [35[. 

The static Domain Wall solution given in [13| has a 
-B(£) oc £, which vanishes at the wall position. This favors 
to impose the 'flat' condition on 4-dimensional fermions, 
as explained on sections !!! A ll and lll A 2l below. However, 
in multi-brane scenarios [l| or in curved branes (23|, a 
different 4-dimensional condition could be imposed, as 
mentioned in sections III A 31 and III A 41 below. 

From now on, the 5-dimensional spinor is decomposed 
into, 



*(z,0 = VM-(z)/ + (0 + lMz)/-(0, 



(4) 



with tp± = P±ip the 4-dimensional chiral spinors and f± 
their profile through the extra dimension. 



1. Flat Massless 4D Fermions 

For a flat brane world volume, the Dirac equations for 
4D chiral massless fermions, ip±, are 



^ (4) V>± = o. 



(5) 



Following the procedure described in one gets 

the equations of motion for the profiles, 

f' ± +A'f±^\V(ct>)e c f±=0, (6) 

where f± are the profiles for ip±. Their solutions are, 



f± oc e 



-A±\ f <H T(<l>)e c 



(7) 



2. Flat Massive 4D Fermions 

If the world volume of the brane is flat, but the 4- 
dimensional fermions are massive, the Dirac equation for 
these are, 



with m the 4-dimensional mass. A similar procedure to 
the one above, yields the profile equations (see fl^|). 

f± + A'f± T W(cf>)e c f ± = ±me- A+c f T , (9) 

since massive fermions mix chiralities, the equations are 
coupled. 

The best way of work them out is changing the variable 
from £ — > £' where £' = J d^er A+c . Then, in terms of 
the new variable one gets 

(10) 

with prime denoting derivative with respect to 

When the equations are decoupled, a Schrodinger-like 
equation governs the profiles (3d| . 

(11) 



f'± + A'f± T XP(4>)e A f ± = ±m/ T , 



[-%> + v£\ 



where 



and 



U± =771 U±, 



f ± i ^ e~ 2A u ± , 



V± m = (XP{4>)e A ) ±d e (\T(cb)e A ) 



3. Covariantly Massless 4D Fermions 



(12) 



There exist branes which are not flat submanifolds em- 
bedded on the spacetime [23|. For these general cases, the 
Dirac equation for 4-dimensional massless chiral fermions 
is given by 

y (4 V± = o, 



(13) 



and the profile equations are, 
1 



f' ± + - A (4A' + 3B') /± T XV{^)e c f ± = 0, (14) 



whose solutions are, 

f± OC e ~i(^+3B)±XfdiV(4,)e c _ 
4- Covariantly Massive 4D Fermions 



(15) 



(8) 



On these curved branes, 4-dimensional massive 
fermions might be considered. Their Dirac equations are, 

)J7 (4) ^ ± = -myb T , (16) 

whilst the profile equations are, 

f± + \ (4A' + 373') /± T XP{4>)e c f ± = ±me c f T . (17) 

Proceeding as before, but with the changes £' = J d^e c 

and f± i — y e i u±, the Schrodinger-like equation gov- 
erning the profiles is, 

(18) 
(19) 



with 



[-<9f, + V^] tt± = m 2 u ± , 



vL = (AP(ci)) 2 ± % {xvm . 
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Remark. Equation ([19]) has the same structure as (TT21) . 
but the lack of the exponential functions make it looks 
simpler. However, non-flat domain- wall solutions are 
much more complicated than flat ones. 



B. Massive 5D Fermions 

Although in theories with extra dimensions it is pre- 
ferred to begin with massless fermions, there exists no 
fundamental reason for doing so. Thus, in case 5- 
dimensional massive fermions are considered, the above 
solutions are valid after the substitution XV ((f)) —> 
XP((f) — M, where M in the 5-dimensional mass. 

From ((7|) and (jT5J) , it is clear that after the substitution 
XP(4>) — > XP((j>) — M, the localization of the positive 
profile, /_)_, is favored in detriment of the negative one, 
/_ . It is possible that a critical region exists where both 
chiralities might be localized. 

In the case of 4-dimensional massive fermions, the ef- 
fect enters through the quantum mechanical potential, 
(TT2"]) and (|19[) . which is a case less intuitive to analyze. 
However, these equations might be solved numerically. 
From their solutions one can argue whether or not the 
model is phenomenologically relevant. This analysis in 
done for the double domain walls in the next section. 



III. DOUBLE DOMAIN WALL 

In the special case of the so called double walls 13], 
whose metric has a line element given by 



ds 2 (g) = e 2A ^ (ri^dx»dx u + d?) 
with a warp factor, 

f 



the scalar field is, 



2s 



ln(l + K) 2s ) 



, ... V6s — 3 , . f 

<p(Q = arctan(a4) 



(20) 



(21) 



(22) 



The quantum mechanical potentials, for 5-dimensional 
massless fermions, and several choices of the polynomial 
V(<f>), are shown in figures[l]|3]and|4] In what follows, the 
parameter a has been set to one and the solid(dashed) 
line indicates the quantum mechanical potential for neg- 
ative(positive) chiralities, as function of the extra coor- 
dinate £. 



A. Massless 5-Dimensional Fermions 

Following the asymptotic behavior of the profiles, 

/±~e^« s («\ (23) 
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FIG. 1: Quantum Mechanical potential of profiles on a DDW 
for linear scalar- fermion coupling. For s £ {1, 3, 5} 

as proposed in 19], it follows that, 



A > 2a 



2s 



TTy/6s — 3 



(24) 



with I the minimum (maximum) power of the polynomial 
P(<f>) if s < 13 (s > 15). 

For the plots presented above, A is set to 



A = 3a 



2s 



TTy/6s — 3 



(25) 



unless it is specified otherwise. 

For linear scalar-fermion coupling, from figure [T] it is 
clear that the positive chirality fermion cannot be local- 
ized for s = 1. However, for higher values of the coupling 
constant [37j , a local well-potential appears around £ = 
(see figured]). The well gets steeper as the Yukawa cou- 
pling increases. 



A very interesting consequence of this is that the ~ 10 6 
factor of the ratio between the Yukawa couplings of elec- 
tron and neutrino could explain why right-handed neu- 
trinos are not seen, while right-handed electrons are. The 
argument is as follows. The mass difference between elec- 
trons and neutrinos, explained by the Higgs mechanism, 
is due to the difference of their Yukawa couplings. As- 
sume that right-handed fermions have the same mass as 
their left-handed partner. Independently of the Yukawa 
coupling, left handed fermions are localized because of 
the well potential V^(£). However, V^(£) for small 
Yukawa coupling (neutrino case) is a barrier potential 
which does not allow localization, whilst for large values 
of the Yukawa coupling (electron case), the appearance 
of a local minimum around the brane position favors the 
localization. 

For quadratic scalar-fermion coupling, the quantum 
mechanical potential is not Z2-symmetric, so it is not 
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FIG. 2: Quantum Mechanical potential of positive chirality 
profile on a Double DW for linear scalar-fermion coupling, 
s—1 and A = 1 (solid-), 2 (dashed-) and 3 (dot-dashed line). 
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FIG. 4: Quantum Mechanical potential of profiles on a DDW 
for cubic scalar-fermion coupling. For s £ {1, 3, 5} 
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FIG. 3: Quantum Mechanical potential of profiles on a DDW 
for quadratic scalar-fermion coupling. For s £ {1, 3, 5} 



physically relevant because the initial problem is clearly 
symmetric under the exchange £ — > — £. The same ar- 
gument can be applied to any even power of the scalar- 
fermion coupling, so they can be excluded from the anal- 
ysis. 

The analysis of the 'spectrum' for generalized odd 
scalar-fermion coupling is quite similar to the linear one, 
except that the "hill" potential, which forbids the local- 
ization of fermions for A < A c , is not present. 

A more detailed analysis of the spectrum of massless 
fermions in these family of domain wall(38j was presented 
recently in [20j | . 



Massive 5-Dimensional Fermions 



On this section a numerical analysis of the moduli 
space is done. Such analysis allows to determine, point- 
by-point on the moduli, whether a tachyon state exists. 
Existence of tachyonic states turns into a phenomenolog- 
ical instability of the theory, which permits to rules out 
such a region of the moduli space, as part of a physical 
model. 

The numerical analysis is done for the action (|3]), with 
A > 1.44, which assures well-behaved asymptotic limits 
for left-handed fermions, as determined in [19J. Both 
chiralities are analyzed independently. 

For solving the profile equations, a Numerov algorithm 
was implemented in Python and Sage [24j . It proceeded 
as follows. First, a lattice was constructed in the mod- 
uli space. The considered lattice covers A £ [1.45,5.00] 
in steps of 0.05, and M £ [0.00, 1.00] in steps of 0.02. 
Secondly, the program finds V^(£' = 0), and from 
those reference start eigen-values (which represents the 
4-dimensional mass square), solves the Schrodinger-like 
equation for the respective "crural" profile. Then, the 
value of the 4-dimensional mass square is swept by in- 
crementing it in 0.2 units, until a range of Am 2 ~ 0.2 
is found, where the ground state mass eigen-value lies. 
Next, the bisection algorithm is applied for improving 
the precision of the 4-dimensional mass-square value. Fi- 
nally, when the 4-dimensional mass-square for the ground 
state is found with an enhanced precision, one might con- 
clude whether it is a tachyonic state, at that point in the 
moduli. 

The scalar-fermion coupling considered below are 
monomial functions of cj), say, V(<f>) = 0*, with t £ 
{1, 3, 5}. Besides, the parameter of the double walls takes 
similar values. 
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TABLE I: Resume of the moduli space analysis. 
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1. Linear Coupling 

When the scalar- fermion coupling is linear (t = 1), 
the value of A should be less than 3.14 ± 0.5 for any 
value of M, in order for the left-handed ground state to 
be physically possible. The right-handed mode does not 
impose any constrain on A. 

Changes in the s-parameter induce changes in X ma x, 
which also varies slowly with M. Additionally, for large 
s, X max increases with M. 

As was mentioned above, there exists an energy gap 
between the left- and right-handed ground states. 

Interestingly, in the simplest case, s = 1, there is no 
possibility of ruling out the tachyonic state for the right- 
handed chirality. This fact somehow points to a non- 
minimalist model for phenomenology. All other values of 
s possess a tachyon free spectrum. 



2. Cubic Coupling 

In the cubic coupling case, V((f) = (f> 3 , the right- 
handed profiles impose no constrain on the moduli space, 
whilst for s = {1,3} the tachyonic state is always present. 

For s = 5 there exists a small window in the mod- 
uli space which allows a tachyon free spectrum, given 
roughly by X max > 1-64. 



3. Fifth Coupling 

When the fifth coupling was considered, besides the 
fact that the right-handed profile imposes no constrain on 
the moduli space, it was not possible to get a tachyonic 
free spectrum for the left-handed state. 



C. Massive 5D Fermions (Revised) 

In odd number of dimensions, it is well known that 
there is a parity anomaly [13, HH ■ In order to avoid this 
anomaly, the fermion mass in odd dimensions should be 
multiplied by a sign function of the scalar field 29] . Thus, 



the action is, 

S = J d 4 x d£ a/IsI* (- y -Msgn(<f>) + XP{4>)) * = 0. 

A numerical analysis similar to the one in the preced- 
ing subsection is done, but considering the action (|26|) 
instead of ©. 



1. Negative Chiralities 

For negative chiralities profiles, once more, the case 
with t = s = 1 has a different behavior than the others. 
Tachyonic states were found all over the moduli space, 
except for M = and any value of A, or for M = 0.02 if 
A > 4.65. 

Any other choice of the parameters M and A gives 
tachyon-free models. It is worth mentioning that choices 
(t, s) G {(1, 3), (1, 5), (3, 1)} show ground states with m 2 
close to zero (m 2 ~ 10~ 6 — 10~ 8 ), although this value is 
far from the numerical error (Am 2 ~ 10 -14 ). 

Due to what seems to be numerical problems in the 
precision, the calculations with t = 5 could not be done. 



2. Positive Chiralities 

On the other hand, positive chirality profiles present 
no tachyonic states, at least for the considered values of 
s, t, M nor A. Thus, the ground state has m 2 > 0, but its 
localization is meta-stable. Besides the fact that meta- 
stable states rise, the question of whether these are long- 
lived enough, has been left aside because the main line of 
this work rests on formal rather than phenomenological 
aspects of the theory [39j]. 

IV. DISCUSSION AND CONCLUSION 

First of all, a new kind of condition on the 4- 
dimensional fermions was considered, 

/ V± = 0, (27) 

for finding the equations for the localized profiles, /±(£)- 
Although only the conditions for localizing modes were 
given, this might be important when multi-brane worlds 
are considered, as well as non-flat branes. 

The developing of a relationship between the 
AdS/CFT correspondence and brane- worlds (see [32j and 
references there in) could add importance to this partic- 
ular sort of conditions, since the correspondence applies 
for a stack of N D-branes. Moreover, in this scenario, the 
symmetry breaking might be interpreted geometrically as 
parallel branes slightly separated in the extra-dimension. 
This last set-up is exactly one kind of scenario where the 
new conditions might be important. 
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Generalizations of the profile equations were found, by 
allowing the scalar-fermion coupling to be non-linear, say 
a polynomial 7 7 (</>) 1 S r \f r . When V(<t>) is a monomial, the 
physically relevant powers are odd [40j, because for even 
powers on <f> the quantum mechanical potential is not 
symmetric around the hyperplane £ = 0. 

Polynomial couplings get complicated as soon as bi- 
nomials are considered, due to the scalar coefficients 
which regulate the form of the quantum mechanical po- 
tential. Focused on the double domain wall family, where 
s £ 2Z+1 labels the members, the condition for localizing 
one of the chiralities is 



A > 2a 



2s 



(28) 



with I the minimum (maximum) power of the polynomial 
P(4>) if s < 13 (s > 15). 

When binomial scalar-fermion coupling are considered, 



V{4>) =e l <f>*< +c 2 4> t >, 



(29) 



it is straightforward to note that both i< and t> should 
be even if the symmetric single brane scenario is exam- 
ined. However, there is a chance that odd powers of <j) 
could be phenomenologically relevant for asymmetric or 
multi-brane solutions. 

Some constraints have been obtain for localizing both 
chiralities of 4-dimensional fermions coming from massive 
or massless 5-dimensional ones. It was found that the 
simplest model, s = t = 1, does not accept localization 
of right chiralities, if the action © is consider. However, 
when analyzing the action (|26l) . both chiralities could 
possibly be localized, and it would provide a natural ex- 
planation of the difference between the left- and right- 
handed neutrino mass, in those models beyond standard 
one where right-handed neutrinos are likely to exist, as 
well as the mass difference between different families of 
fermions. For non-linear scalar-fermion the tachyon-less 
constrain is made more restrictive as the non-linearity 
increases. 

When a 5-dimensional massive fermion is considered, 
a slight correlation between M and A was found, in the 
studied region (0.00 < M < 1.00 and 1.45 < A < 5.00). 
The considered region might be extended, but covering a 
larger moduli space goes beyond the scope of this work. 

One might wonder that even when no additional gauge 
structure is given, left- and right- handed chiralities be- 
have differently. An important question to be answered 
is whether charged particles could have similar behavior 
for both chiralities whilst neutrals do not, because par- 
ticles other than neutrinos have the same mass for both 
chiralities. In order to follow this line a gauge field, say 
Abelian one, must be taken into account. This proposal 
could be part of future research. 

Similar analysis to the one done in this work might 
be done for the asymmetric D-branes in [l3|-[l5j], some 
advances towards this are already being done by the au- 
thors. 



In [3l|, the behavior of fermions on (anti-)kink solu- 
tions under discrete symmetries is studied. These could 
be used to complete a theoretical analysis on fermion lo- 
calization. 
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Appendix A: Conventions and Notations 

Through the manuscript, the metrics have signature 
mostly positive. Since the tetrads formalism is used ex- 
tensively, distinction between flat and curved coordinates 
is made by Latin and Greek indices. Moreover, hated in- 
dices run over the whole spacetime whilst unhated ones 
run over a hypersurface restriction, i.e., on the coordi- 
nates parallel to the topological defect. 

The gamma matrices are defined in the tangent space, 
and they satisfy the Clifford algebra, 



{ 7 a ,7 S } = 277 aS l. 



(Al) 



In even dimensions one can define the chirality matrix 
7*, satisfying the properties 

• {7 S ,7*} =0. 

• (7*) 2 = 1- 

From this, the projector operators, 

1-7* 



(A2) 



are both non-trivial. 

In any dimension one may define a set of generators of 
the Lorentz algebra, constructed with the elements of the 
Clifford algebra (jAll) . These generators of the Lorentz 
algebra are, 



7 ,7 



(A3) 



When the Dirac-Feynman slash notation is used it 
must be interpreted as, 



9= 



(A4) 



by. 



The covariant derivative for gauge theories is defined 



V/j = dp, - igAp, - iVtp, 



(A5) 
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with the gravitational connection, which is related to 
the Christoffel connection for integer spin fields, and with 
the spin connection for semi-integer spin fields. Clearly 
the Dirac-Feynman slash notation can be used with the 
covariant derivative, 

7= £ Q A 7 3 V A = El 7 a - igAt - *0 A ) . (A6) 
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